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Abstract 

Large-scale universal quantum computing requires the implementation of quantum error cor- 
rection (QEC). While the implementation of QEC has already been demonstrated for quantum 
memories, reliable quantum computing requires also the application of nontrivial logical gate op- 
erations to the encoded qubits. Here, we present examples of such operations by implementing, in 
addition to the identity operation, the NOT and the Hadamard gate to a logical qubit encoded in 
a five qubit system that allows correction of arbitrary single qubit errors. We perform quantum 
process tomography of the encoded gate operations, demonstrate the successful correction of all 
possible single qubit errors and measure the fidelity of the encoded logical gate operations. 

PACS numbers: 03.67.Pp,03.67.Lx 
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FIG. 1: Outline of the quantum algorithm. 

Introduction. -Quantum computers can solve certain problems exponentially faster than clas- 
sical computers jl, 2|. An essential precondition for realizing this potential is the preservation 
of the coherence between quantum states. This requirement makes the implementation of 
quantum computing much more challenging than for classical devices. Reliable quantum 
computing ^| is possible, in principle, provided quantum error correction (QEC) schemes 
can be implemented with a fidelity above a certain threshold [l[ \^L\- Every QEC code has 
an overhead in terms of gate operations and additional (ancilla) qubits. The protection of a 
single qubit against arbitrary single-qubit errors requires at least five physical qubits js, 9]. 
Over the last few years, a series of experiments were performed that demonstrated that 
QEC is indeed capable of protecting quantum states against generated (artificial) errors 
or decoherence induced by the environment (see, e.g. 10Ml4|). However, the realization of 
reliable quantum computing requires more than the preservation of information: it must also 
be possible to process the encoded information by applying logical gate operations to the 
protected qubits. Here, we present an experimental demonstration of such gate operations 
on an encoded qubit. For encoding, we use a five bit QEC code jl], [s| that allows correction 
of arbitrary single-qubit errors (a so-called perfect QEC code) and demonstrate successful 
gate implementation and error correction. The result is a fault-tolerant implementation of 
the corresponding gate operations. Fig. [1] summarizes the scheme: it starts with encoding 
the input state of the first qubit into five physical qubits. To the resulting state, we apply 
one of three single-qubit gates - the identity, NOT or the Hadamard gate. As the third step, 
we apply another operation, which is either the identity (corresponding to no error) or one 
of the fifteen possible single-qubit errors. In the fourth step, the information is decoded, i.e. 
the output is extracted into the state of the first physical qubit. In the fifth and final step, 
possible errors are detected and corrected. 

For the physical qubits, we use a system of five nuclear spins. The molecule containing 
the spins is dissolved in an anisotropic solvent. The resulting magnetic dipole couplings 
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between the nuclear spins are significantly stronger than the more frequently used scalar 
couplings and therefore result in a speedup of the gate operations by approximately an 
order of magnitude. The complete operation can therefore be completed within a time 
period significantly shorter than the coherence time of the system. 

Five-qubit error correcting code. - Fig. [2] shows the quantum circuit for the five-qubit QEC 
code [8|. Qubit 1 is the register bit carrying the input state \ip) = a\0) + j3\l), where a and 
(3 are arbitrary complex numbers with |a| 2 + |/3| 2 = l. The other four qubits are initialized 
in the state jO)® 4 . The unitary operation U en implemented by the circuit encodes the state 
| <f) into a logical state as 

U en (a\0) + /3|1»|0000> = a\0) L + (3\1) L , (1) 



where 



\o) L = 4=(|ooooo) - |10111) - |01011) + |11100) 



|10010) + |00101) + |11001) + |01110)) (2) 



\1) L = -^=(|11111) - |01000) + |10100) - |00011) 



+ |01101) + |11010) - |00110) - |10001)) (3) 

are the computational basis states of the logical qubit. The decoding operation Ude is the 
inverse operation of U en , i.e., Ude = U\ n . 

The five-qubit QEC code can detect and correct arbitrary single-qubit errors. The possible 
single-qubit errors for a five-qubit system can be written as bit flip errors Bk, phase flip 
errors Sk and combined bit- and phase flip errors BSA;, where k = 1 ... 5 indicates the 
affected qubit. These fifteen errors, together with the identity operation E define the possible 
outcomes if only single-qubit errors occur. The 2 4 possible states of the four syndrome qubits 
can distinguish between these 16 different outcomes. This is used by the error correction 
step, which is a unitary operation on the first qubit, controlled by all four syndrome qubits. 
Without the encoded gate o per ations, this code was implemented previously in a system of 
weakly coupled spin qubits [11 ]. 

Experimental protocol. -For the experimental implementation, we used the two fluorine and 
three proton spins of the molecule of l,2-difluoro-4-iodobenzene, whose structure is shown 
in Fig. [3j The molecule was dissolved in the liquid-crystal solvent ZLI-1132 to retain the 
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FIG. 2: Network representation of the encoding operation. The input state \(p) = a\0) + /3\1) is 
initially encoded in qubit 1. Qubits 2-5 are the syndrome bits. H denotes a Hadamard transform, 
and Z the tt phase gate. The box or (J) with a connected vertical line indicate controlled gates, 
where the filled or empty circle marks the control qubit. The operations are conditional on the 
control qubit being in state |1) (filled circle) or |0) (empty circle), respectively. The output of the 
circuit is the five-qubit state ol\Q)l + 
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FIG. 3: Parameters of the spin qubits in the molecule l,2-difluoro-4-iodobenzene. The inset shows 
the structure, where the five qubits 1-5 are spins Fl, F2, HI, H2, and H3, respectively. The 
chemical shifts with respect to the transmitter frequencies of the proton and fluorine spins are 
shown as the diagonal terms and the dipolar couplings between spins are shown as the off-diagonal 
terms in units of Hz. The effective relaxation times T| are determined by fitting the peaks in the 
spectra. 

dipolar couplings between the spins. We denote the two fluorine spins Fl and F2 and the 
three protons HI, H2, and H3 as qubits 1-5. Data were taken with a Bruker Avance II 
500 MHz spectrometer. The relevant Hamiltonian of the dipolar coupled spin system is, in 
frequency units, H = J2i + Si<j ^-tji wnere HZ — ^ViZi is the Zeeman Hamiltonian and 
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TLfj = irDijZiZj describes the dipolar coupling Hamiltonian for the heteronuclear spins % and 
j, or TLfj = (irDij /2)(2ZiZj — XiXj — YiYj) for the homonuclear case. The indices i and j run 
over all five spins, Xj, Yi, Z^ denote the Pauli matrices, the chemical shift of spin i, and 
Dij the coupling constant. Compared to the dipolar couplings, the scalar couplings between 
the spins are negligibly small or can be merged in the dipolar couplings 15) . To determine 
the numerical values of the Hamiltonian parameters, we measured different experimental 



spectra, using multiple quantum NMR 



16| and heteronuclear decoupling. These spectra 



were used in a fitting process which yielded the parameters listed in Fig. [3j 
We prepared the initial states X0000 or F0000, with = |0)(0| by the circuit in Ref. (lj 
The required unitary operations were implemented by two optimized shaped pulses designed 



with a GRAPE algorithm 
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19] . The experimental spectrum of the state X0000 is 



It contains a single main peak, which is 



illustrated in the Supplementary Material [15 
the signature of a pseudopure state (PPS) (l7|. By comparing with a reference spectrum 
obtained from the thermal state, we estimate that the polarization of this state is ~ 0.72 
of the maximum polarization that could result from an ideal preparation. For the following 
experiments, we normalized the spectra to this one, so the fidelities of the QEC protocol 



do not include losses during the PPS preparation [15]. The Z0000 state was prepared by 
applying an additional tt/2 rotation to the X0000 state. After the encoding step, we applied 
one of the following gate operations to the logical qubit: the identity E, the NOT gate (up 
to a known phase), or the Hadamard gate. 

The NOT gate Nl for the logical qubit is relatively simple to implement for this code, since 
it is transversal, i.e. it can be written as 



N L = Ry(nf 



(4) 



where R y (n) = e lnY / 2 . This relation can be verified by comparing the definitions of the 
logical states in Eqs. (T5H3]). The elements of N L in the logical basis are represented as 
(0l|A l |0 l ) = (1 L |JV £ |1 L > = 0, (1 L \N L \0 L ) = i, and (0 L \N L \1 L ) = -%. 

In contrast to the NOT gate, the Hadamard gate is not transversal in this code [20|. We 
therefore have to design an operation that implements this gate in the 2 5 dimensional Hilbert 
space of the five-qubit system. In this space, the Hadamard gate should generate a tt rotation 
around the (1,0,1) axis of the two-dimensional subspace spanned by the states \0)l and 
and an identity operation on the other 30 states. The corresponding unitary operator thus 
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has the matrix representation 





1 
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V2 


V2 



Hr 



(5) 



V i/ 

where zero elements are not shown. This unitary operation can again be implemented by 
an optimized shaped pulse, designed in the same way as the other pulses described above. 
After the unitary gate operation, we also applied Bk, BS/c, or Sk errors to the individual 
qubits, using single-qubit n rotations around the x, y, or z-axis. 

For a quantitative evaluation of the algorithm's performance, we used quantum process 
tomography (QPT) 2l| of the complete algorithm represented in Fig. [TJ The process 
can be completely characterized by its x matrix which maps an arbitrary input state 
p in into the output state p out = ^ w Xwe fc /o in eJ. Here the operators e k ,i G {E,X,-iY,Z} 
denote the basis set for describing the process, and the indices k, I — 1, 4 run over the 
elements of the basis set. The measurement of \ requires the preparation of four input 
states pi n = E, X, Y, and Z. For each p in , we determined the output state p out by quantum 
state tomography. Since the unit operator E is always time- independent, the corresponding 
input state is omitted, assuming the output state is E. The quantum state tomography is 
performed by measuring in one experiment the transverse magnetization and in a second 
experiment applying a 7r/2 readout pulse and then measuring the transverse magnetization. 
The measured FIDs were Fourier transformed and the resulting spectra were fitted to the 
theoretical spectra by adjusting as a single parameter the overall amplitude. Having the 
output states for the four input states, we determined the the x matrix using the established 
strategy For each process, we quantified the performance by comparing the experimental 
iXexp) and theoretical (xth) X matrices via the fidelity [22] 



F X = \ Tl (XexpXth)\N^(XexpXexp)Tl(XthXth)- 



(6) 



Experimental results. -We first checked the encoding, encoded gate and decoding operations 
with a simplified experimental scheme. Compared to the full scheme shown in Fig. [TJ 
we omitted the the error and error correction operations. Fig. H] shows the experimental 
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FIG. 4: Bar plots of the real parts of the x matrices for the identity, NOT and Hadamard gates in 
logical states. The RMS values of the imaginary parts are 0.035, 0.0069, and 0.017, respectively. 

results for the process matrices \ as bar plots. The corresponding matrices for the ideal 
operations have xn — 1 f° r the identity operation, ^33 = 1 for the NOT operation, and 
X22 = X44 = X24 = X42 = 0.5 for the Hadamard gate. All other elements should vanish. 
We find a good qualitative agreement between the theoretical and experimental values, with 
fidelities of 0.979, 0.983, and 0.956, for the identity, NOT and Hadamard gates, respectively. 
We next tested the full implementation, including errors and error correction. Fig. [5] il- 
lustrates the results for three out of 48 experiments. In the first example (trace b), we 
initialized the system to the Y state and applied an E gate and SI error operation, followed 
by the decoding and error correction steps. Trace c) shows the corresponding results for the 
X initial state, NOT gate and BS4 error and d) for Z initial state, H gate and B5 error. 
The insets show enlarged partial spectra containing the main signal components, with the 
experimental spectra represented by dashed lines, the ideal spectra as full lines. Experimen- 
tal and theoretical curves were both normalized to the spectra of the initial PPS. Figure [6] 
shows the measured fidelities for each of the 48 different experiments as a bar plot. The 
solid horizontal line shows the average fidelity for each type of gate, averaged over the 16 
different error conditions. 

To assess the usefulness of the scheme, we compare the achieved fidelities to an idealized 
experiment where we do not use QEC, but the same 16 error conditions can occur, with equal 
probabilities. In this case, the three single qubit errors acting on the first qubit result in zero 
fidelity, while the other 13 error conditions (the identity and the single-qubit errors on the 
ancilla qubits) result in fidelities of one. Averaged over these 16 reference experiments, we 
would thus expect an average fidelity of 13/16 = 0.8125. This value is shown in the figure 
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FIG. 5: (Color online.) Spectra of the fluorine spins obtained for different versions of the experi- 
ment, a) Reference spectrum obtained by applying a tt/2 pulse to the thermal equilibrium state, 
b) Spectrum of the output state, starting with the initial state Y, encoding, applying the identity 
operation, SI error and error correction, c) Same as b), but pi n = X, NOT gate and BS4 error, d) 
Same as b), but pi n = Z, Hadamard gate and B5 error. The dashed and solid curves indicate the 
experimental and simulated spectra, which are normalized to the initial pseudopure states. The 
lines with the main signal contributions are shown enlarged as insets. 

as the horizontal dashed line. The average fidelity for the three different gates exceeds 
this reference value by 0.0837, 0.0528 and 0.0196, for the identity, NOT and Hadamard 
gates, respectively. This shows that the performance of the QEC scheme is high enough to 
compensate the additional errors affecting the syndrome qubits as well as the errors due to 
experimental imperfections of the encoding, decoding and error correction steps. 
Discussion. -The fidelities achieved in these demonstration experiments are still lower than 
the threshold fidelities required for scalable quantum computation. We have identified four 
main causes for the observed reduction of fidelity: (i) Stability and homogeneity of the 
magnetic and radio frequency (r.f.) fields were less than ideal, since the lack of deuterium 
in the sample did not allow operation of the lock system of the spectrometer, (ii) Finite 
accuracy of the Hamiltonian parameters, (iii) Finite coherence time: the duration of the 
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FIG. 6: Bar plots showing the experimentally determined fidelities for the identity, NOT and 
Hadamard gates applied to the encoded qubit after applying error operations and error correction. 
The average fidelity is shown as the solid line. Compared to the average fidelity without error 
correction, shown as the dashed line, the encoding - decoding - error correction scheme improves 
the fidelity. 

experiment is in the range of 26-64 ms (excluding the preparation of the initial PPS), which 
is comparable to the measured T|. (iv) Deviations between the calculated pulse shapes and 
those acting on the spins, which are caused mostly by nonlinearities of the spectrometer 
hardware. Our efforts to suppress these imperfections concentrated on 1) using short gate 
durations, 2) designing the pulses robust to variations of the r.f. strength and insensitive to 
frequency offsets and 3) 'fixing' the pulses by measuring the actual amplitudes of the control 
gates with a pick up coil and adjusting the amplitudes to minimize the difference between 



theoretical and experimental values 23]. 



For a better understanding of the sources of errors, we performed a quantitative analysis 
of one of the experiments: the Hadamard gate combined with a phase flip error at qubit 5 
and the corresponding error correction 15j. Combining experimental results and simulation, 
we estimate that the field inhomogeneity and the imperfect implementation of the pulses 
contribute ~ 0.09 to the loss of fidelity. Additionally the limit of T 2 and the imprecision in 



characterizing the Hamiltonian, contribute ~ 0.04 and 0.03 to the loss of fidelity, respectively. 
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The simulation with results in a fidelity below the experimental value, indicating that the 
coherence time in the experiment is longer than T 2 *. This is expected if the gate operations 
refocus some of the inhomogeneity of the system. 

Conclusion. -The results presented here are a first demonstration of one of the most impor- 
tant preconditions for reliable quantum computation: the implementation of logical gate 
operations on encoded qubits. While QEC has been demonstrated on quantum memories, 
its combination with processing of quantum information is an important milestone towards 
the implementation of reliable quantum computing. The present demonstration used three 
single-qubit gate operations acting on a logical qubit encoded in a perfect five-qubit QEC 
code. In the earlier implementation of this code ll|, without the gate operations, the ex- 
perimental duration was > 300 ms. Compared to that, we have reduced the duration of the 
experiment by approximately one order of magnitude by using a dipolar coupled system, 
whose interactions are stronger than the scalar couplings used before. The reduction of the 
gate operation times opened the possibility to implement in addition the logical gate oper- 
ations. With better homogeneity and more precise Hamiltonian parameters, it should be 
possible to improve the experimental fidelity. The recent theoretical progress in measuring 



Hamiltonians with dipolar couplings 



24[ and optimal algorithms in pulse finding 25] should 



be helpful for this purpose. This will allow us to control larger systems, encoding multiple 
qubits and implementing multi-qubit gate operations, such as CNOT in the encoded qubits. 
If the fidelity can be improved sufficiently to reach the error threshold, the combination of 
QEC with logical gate operations will pave the way to reliable quantum computation. 
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